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ABSTRACT 


The  detailed  development  of  the  analyses  which  are  the  bases  for  the  two- 
stream  axisymmetric  base-pressure  programs  (TSABPP-1,2)  is  presented. 

For  the  "corresponding"  inviscid  flow  field,  these  analyses  include  the  Prandtl- 
Meyer  expansion,  the  oblique  shockwave,  the  slip  line,  and  the  method  of 
characteristics.  For  the  turbulent-mixing  component,  two-dimensional  constant- 
pressure  mixing  and  its  application  to  axisymmetric  flows  are  considered. 
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i.  INTRODUCTION 


The  purpose  of  thi  ■.  report  is  to  present  the  detailed  development  of 
the  principal  relationships  ai  solution  techniques  which  are  the  bases  for  the 
two-stream  axisymmetric  base-pressure  programs  (TSABPP-1,  2)  .  While  the 
intent  here  is  to  inform  the  interested  program  user  about  these  specific  pro¬ 
grams,. it  should  also  be  noted  that  various  subsections  of  these  programs  are 
more  general,  and  as  such,  can  serve  as  basic  building  blocks  in  the  analysis 
of  other  problems.  Although  an  attempt  has  been  made  to  make  this  report 
complete  in  itself,  it  is  intended  and  should  be  considered  as  supplemental  to 
the  technical  reports  [1,  2]  documenting  the  computer  programs. 

A  uniform  notation  has  been  maintained  through  the  reports  and  the  com¬ 
puter  programs;  since  comprehensive  lists  of  symbols  are  included  in  the  above- 
referenced  technical  reports,  they  are  not  repeated  here. 


2.  DISCUSSION 

In  the  development  of  the  base-pressure  programs,  the  analyses  and 
the  resulting  programs  were  subdivided  into  three  major  parts.  The  first  part 
was  the  calculation  of  the  "corresponding"  inviscid  flow  field;  the  second  part 
was  the  calculation  of  the  mixing  component;  and  the  third  part  was  the  overall 
organization  of  the  former  calculation  sequences  into  a  master  program  for 
determining  by  iteration  the  base-pressure  and  base -temperature  solutions. 

This  note  is  not  concerned  with  the  latter  part,  since  it  is  discussed  in  detail 
in  [1],  but  rather  with  the  former  parts. 

The  inviscid  flow-field  analysis  consists  of  the  calculation  of  the  < 

supersonic  external  (free  stream)  and  the  internal  (nozzle)  flows  b}'  the  method 
of  characteristics.  The  external  flow  is  over  a  solid  afterbody  followed  by  a 
constant-pressure  boundary  region.  The  internal  flow  is  assumed  to  be  from 
either  an  ideal  conical  or  uniform  flow  nozzle  discharging  into  a  region  at  con¬ 
stant  pressure.  The  supersonic  external  and  internal  flows  interact  at  their 
impingement  point,  if  it  exists,  to  form  an  oblique  shock  system  and  a  slip  line. 

The  component  flow-field  analyses  are  discussed  in  Appendixes  A  and  B. 

The  mass  entrainment  and  energy  transport  rates,  due  to  the  turbulent 
mixing  regions  formed  between  the  fluid  in  the  wake  region  and  the  internal  and 
external  streams,  are  estimated  by  locally  superimposing  the  two-dimensional 
constant-pressure  turbulent-mixing  model  at  the  impingement  point  of  the  "cor¬ 
responding"  inviscid  streams.  The  mixing  analysis  and  the  "corresponding" 
inviscid  flow-field  analyses  are  then  related  by  the  recompression  criteria  and 
the  requirement  of  conservation  of  mass  and  energy  in  the  wake  region.  These 
analyses  are  discussed  in  Appendixes  C  and  D. 
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Appendix  A 

MISCELLANEOUS  GAS  DYNAMICS  FUNCTIONS 


1.  Prandtl-Meyer  Expansion  (EMSPM) 

The  flow  properties  before  and  after  a  centered  Prandtl-Meyer 
expansion  (Figure  A  — 1)  are  related  through  the  following  expression  (13  ] : 


w«M»,  M* 

\  2  max 


;)  =  “  (Mr 


*,  M*  +  6 
1  max  / 


( A—l) 


where  the  turning  angle,  <5,  is  defined  as 


6  =  T(02  -  00 


for  left- or  right- running  waves.  The  u>  (M::,  M*  \  funetion  is  defined 

y  max J 

)  =  Ivl*  tan-1  ,  — M>f  1 —  *  , 

\  niaxy  max  m*2  -  M * 2 

max 


-i  .,...2  M*  -  1 

tan  H"  - - - - 

1  max  ... .  2 


M*  -  M* 
max 


where  ^corresponding  to  M  —  M*  — 


“L=  l  (T  +  1)/<Y  *  1)  1' 


(A-4) 


The  usual  problem  to  be  solved  is  with  values  of  [Mf,6,y]  specified,  to 
find  from  equations  (A-l)  and  (A-3)  the  solution  value  of  M2‘.  The  numerical 
solution  of  these  equations  can  be  easily  and  quickly  accomplished  as  follows. 

For  expansions,  the  value  of  |  must  be  in  the  range 

y  2  max/  b 

«(M*  M*  \  w(M*,  M*  ^  s  ^  [m-  -  ll  (A -5) 

y  1  max/  y  2  max  j  2  max  J 

and  correspondingly,  the  solution  value  of  M#  must  be  initially  in  the  range 


M*  s  M* 

1  2  max 


(A-G) 


P*M  EXPANSION 


M2  (M|) 


FIGURE  A-l.  PRANDTL- MEYER  EXPANSION  NOTATION 


For  small  turning  angles,  oblique  shock  waves  can  be  approximated  by 
'reversible  compressions";  the  initial  solution  range  for  those  cases  is: 


0  w(M*  M*  \  -  w(M*  M*  \ 
y  2  maxi  y  1  maxi 


(A-7) 


1  -■  s  Mf  . 


(A-8) 


For  either  an  expansion  or  an  admissible  "reversible  compression,"  rapid 
convergence  to  the  solution  value  of  Mjf  is  achieved  by  a  process  of  interval 
halving  and,  at  the  same  time,  successive  reduction  of  the  possible  solution 
interval.  The  solution  is  always  bounded  on  the  negative  or  positive  side  relative 
to  the  sign  of  the  difference 


If  d ^ n_)  <  (),  M^n  ^  is  replaced  by  and  M^a  ^  --  M*  ^  .  On  the  other 

hand,  if  d  ^  ^  >  0,  M*Jn  ^  is  replaced  by  and  Mf,' n  ^  This 

P  N  N 

process  is  continued  until  the  solution  is  isolated  within  a  calculated  interval; 
i.e.  ,  the  sign  of  has  changed  at  least  once.  Convergence  to  the  final 
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solution  is  then  achieved  by  interpolation  and  a  continuing  reduction  of  the  solu¬ 
tion  interval  as  before;  the  solution  procedure  is  illustrated  for  an  expansion  in 
Figure  A-2. 


FIGURE  A-2.  SOLUTION  PROCEDURE  FOR  A  PRANDTL- MEYEIi 
EXPANSION 


The  solution  criteria  are  specified  as 


d(UVw(M*  M* 


max 


< 


t 


and/  or 


^Typically  for  q  = 


tz  =  10  4  <  n  <  G  iterations  are  required. 


o 


SHOCK 


FIGURE  A-3.  OBLIQUE  SHOCK-WAVE  NOTATION 


oblique  shock-wave  angle  (a)  ,  the  approach  Mach  number  (Mt) ,  the  turning 
angle  ( 6) ,  and  the  specific  heat  ratio  of  the  gas  (y)  is  given  [4]  as: 

{sura)3  +  b(sin2a)2  +  c(sin2  o)  +  d  =  0  (A-9) 


where 


b  = 


y  sin2  fi 
M? 


(A-10) 


^Typically  for 


10-1,  4  <  n  <  6  itei’ations  are  required. 


I! 


and 

y 

[ -B/ 2 j  S  [-A3/27 ]/2  .  (A- 18) 

If  equations  (A-16)  and  (A-17)  are  satisfied,  the  three  roots,  Y(I),  will  all  be 
real.  The  smallest  root  corresponds  to  a  physically  impossible  process;  the 
next  larger  root  corresponds  to  a  weak  oblique  shock;  and  the  largest  root 
corresponds  to  a  strong  oblique  shock. 

After  the  weak-shock  solution  value  of  (sin2  u)  has  been  four"1  all  other 
pertinent  dimensionless  ratios  for  an  oblique  shock  wave  can  be  found 
[4,  equations  128  through  149],  Specifically,  the  static  pressure  ratio  across 
the  shock  is  found  from: 
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(A-19) 


P2_  2y  Mi  ( sin2  a  )  -  (y  -  1) 
Pi  (Y  +  1) 


3.  Slip-Line  Analysis  (SLIP) 

For  an  initial  slip  line  solution  range,  assuming  it  exists,  it  is  known 
that  (Figures  A-4  and  A-5) 

fl2  <  0S  <  ^  (A-20) 

since 

I'Ol 

;  ([  6j  -  6i  “  (Oi  -  0a)  ,  (A-21) 

(°2  -°s)| 

subtracting  equations  (A-21),  the  result  is: 


(6i  -  £j)  =  (01  -  0j)  ■  constant  across  the  shock  system.  (A-22) 
Let  (/^and  (inbound  the  solution  on  the  right  and  left,  respectively.  Initially, 

iS  cS 


•u-  »il 


since  0  <  6  <  01  . 
u  S  X 


(A-23) 


The  solution  range  can  possibly  be  narrowed  for  the  P~0  characteristics  (weak 
solution)  as  shown  in  Figure  A-5.  Following  this  figure,  let 


M  (01  "  ^Im) 


!2M  V  2  +  62m)) 


0 

lM 


(A-24) 


where  5  ,  are  the  maximum  turning  angles  for  the  given  y.  and  M*.  Now  for 
iM  1  i 

the  solution  range, 


a)  “"iM’V  set02s~  “lM 


(A-25 ) 


b)  "m 


a 


FIGURE  A -4.  SLIP  LINE  FLOW-FIELD  CONFIGURATION 


Therefore,  depending  on  the  magnitudes  of  0lf  02.  ^2M*  so*u“ 

tion  range  is  defined  by 


(1) 


°1M  >  V  °2s  ~  °1M 


or 


<  0o, 


■<l)  -t  o 
2s  2 


r°l8  =  °r  ei<e2M 


<  u  < 
s 


°ls  °2M’  °1>  °2M 


(A-26) 


For  the  possible  solution  range  <  0  <  0^\  does  a  solution  exist? 

ZS  S  IS 


If 


pi:)<  P.  M  and  pf  1^<  P 
2a  1M  Is  2M 


]■ 


then  a  solution  must  exist  in  the  range 


i/1,  <  0  <  Cases  1  and  2,  Figure  A-6,  illustrate  typical  slip-line  solu- 

2s  s  Is 


tions  and  the  determination  of  the  initial  solution  range. 


FIGURE  A-6.  SOLUTION  EXAMPLES 

If  P^1'  >  P„ „  and/or  p!1'*  >  P„„_,  no  weak  solution  exists  for  the  given 
2s  1M  Is  2M 


conditions  (Case  1,  Figure  A- 7) . 

If  01M  >  ®2M  Nosolution  (Case  2»  Figure  A-7) .  ; 

If  >  0-,j,  -“Possible  solution;  needs  further  analysis  (Cases  1,  2 
in  Figure  A-6  and  Case  1  in  Figure  A-7) . 


FIGURE  A-7.  NO-SOLUTION  EXAMPLES 


AfJ?er  the  existence  of  a  slip-line  solution  and  the  initial  solution  interval, 

,  have  been  established,  convergence  to  the  solution  value  of  the 

slig-Bxie  angle  9  t  can  be  rapidly  achieved  by  a  process  of  interval  halving  and 
s 


successive  reduction  of  the  solution  interval.  The  solution  interval  is  reduced 
according  to  the  sign  of  the  normalized  difference  in  the  pressure  ratio 
(Figure  A-5) , 


,(n) 


.[p;v-pf»]/[o.5(p<",  +  pf))] 


(A-27) 


If  d^n*<0,  then  ^  is  replaced  by  9  ^  and  0,^  =  9^1  1’.  Orifd^>0, 


Is 


2s 


2s 


0  ^  is  replaced  by  0  ,n^  and  =  9 \n  ^ .  This  process  is  continued  until 

*8  S  IS  iS 


it 


the  solution  criterion 

I  ,(n)| 

d  <  e 

is  satisfied.  Typically  for  e  =  10“5,  convergence  is  achieved  for  n 


Appendix  B 

METHOD  OF  CHARACTERISTICS 


This  appendix  is  based  in  part  upon  a  previous  report  [6] . 


1.  Basic  Equations 


For  steady  irrotational  axisymmetric  flow  [3],  the  differential 
equation  for  the  complete  velocity  potential  in  terms  of  the  cylindrical  coor¬ 
dinates  (X,  R)  is 


where 


U 


u2\  32^  uv  32$  /  v2\ 82$  _ 

c2/  ax2  ~  2  c2  axaR  +  v  c2/  sr2  +  r  ~ 


(B-l) 


c2 


(u2+  v2);  V2  =  (u2  +  v2) 


(B-2) 


and 


u 


8$  _  8$ 
ax  '  v  ~  8R  ■ 


The  condition  that  the  derivatives  of: 


3$ 

*3X 


and 


d<P 

dR 


<B-3) 


may  be  discontinuous  along  curves  on  the  solution  surface  to  equation  ( B— 1) 
/  a2<|>  82$  32$ 

I  which  implies  that  — — j  ,  — y ,  and  —  ■ - ■■  are  indeterminate  along  such 

\  a  A  .  oAoJrt 

curves  —  the  characteristic  curves!  yields  the  following: 

a)  The  physical  characteristics: 


(it),  u  -  <"  * al 


(B-4) 
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b)  The  hodograph  characteristics: 


(do] 


I.  II 


cot  o/  sin  0  sin  a 

s~0-a] 


(dR) 


I,  II 


R 


=  0 


(B-5) 


The  Family  I  (right- running)  and  Family  II  (left- running)  characteristics  and 
the  applicable  notation  are  shown  in  Figure  B-l. 


R 


Introduction  of  the  velocity  of  sound  at  sonic  conditions  as  a  reference 


and 


C* 


2yHt 

y  +  1 


(B-6) 


M*  =  V/C* 


(B-7) 


into  equation  (B-5) ,  yields 


(dM*)  rr  ,  .  „  . 

.  ,  I,  II  1  sin  0  sin  a 

I,  II  ±  ( M* tan  a)  T  R  sin  (G  =f  a)  (dR'l,II  "  0 


(B-8) 
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The  problem  of  solving  equation  (B-l)  subject  to  the  imposed  boundary 
conditions  then  becomes  equivalent  to  solving  the  set  of  simultaneous  ordinary 
differential  equations  (B-4)  and  (B-8)  under  the  same  conditions.  Various  "unit 
processes"  encountered  in  the  numerical  solution  of  equations  (B-4)  and  (B-8) 
will  now  be  discussed. 


2.  Field  Points  (FPS) 

For  two  known  points  (1)  and  (2) ,  Figure  B-2,  on  the  Family  1 
and  II  characteristics,  respectively,  the  location  and  flow  variables  at  the 
unknown  point  (3)  at  the  intersection  of  these  characteristics  can  be  determined 
by  use  of  the  characteristics  relationships,  equations  (B-4)  and  (B-8). 


FIGURE  B-2.  GENERAL  FIELD  POINT 


To  a  first  approximation,  the  coordinates  of  point  (3)  are  given  in  finite 
difference  form  as 

_  [(R;  -  Ri)  +  Xt  tan  (0  -  a)a  -  X2  tan  (6  +  a )23j^ 

3  ~  [tan  (0  -  a)  13  -  tan  (0  +  a)  23] 

^ The  notation  (  )  indicates  that  average  values  between  the  points  (ii 

1«J 

and  (3)  are  to  be  used. 
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where  [o  +  ( - 1 ) 1  a  ]  is  defined  as  the  average  value  of  the  quantity  between  the 
1 

points  (i  =  1,2)  and  (3).  C or x-es ponding ly,  the  flow  variables  at  (3)  are  given 
by: 


_  I  i513^1  +  P23,J2+  (Ml  -  Mj>  )  +  Q13(R3  “  K[)  -  Q23(R3  “  R2)l 

lPi3  +  P23  ] 


(B-ll) 


a-Ij  =  [M?  -  p13(ti3  -  Oi)  +  Qi3(R3  -  Ri)J  , 


(B-12) 


where  Q  „  and  P  (i  =  1  or  2)  are  coefficients  based  on.  the  average  values 
id 

between  the  poirn  (i)  and  (3). 


These  coefficients  are  defined  as: 


Pi3=  (M:;:taaa)i3 


(B-13) 


n  u  sin  0  sin  a 

1  R  sin  [0  +  (-1)  a] 


(B-14) 


The  values  of  [X,  R,  M!\  0]3  are  determined  initially  by  assuming  that 
the  flow  variables  at  point  (3)  are  simply  the  average  of  those  at  points  (1)  and 
(2) ;  hence, 

03  =  \  (Oi  +  02) 


RI3  -  2  ( M;  +  M2) 


Then  by  using  equations  (B-9)  through  (B-12)  in  sequential  order,  a  first  approxi 

mation  to  the  values  [X,  R.  M:;‘,  0)3*^  can  be  determined.  These  values  are 
then  used  to  detei’inine  the  avei’age  quantities  and  subsequently  the  next  approxi- 


malion  [X,  R,  Ivl*,  0]  y  .  if  this  procedure  of  successive  approximations  is 
repeated,  values  for  the  variables  at  point  (3)  for  two  successive  approxima¬ 
tions  (n  -  1)  and  (n)  will  be  obtained;  i.e., 


[X,  R,  M*.  0)3(n'1}  and  (X,  R,  AP.w  Js(n)  . 

If  the  problem  areas  discussed  below  are  not  encountered,  the  values  calculated 
by  tliis  procedure  stabilize  rapidly,  and  the  iteration  is  terminated  when 


r  (n)  (n-1)" 

In) 

o 3  -  ^3 

lu  3 

(B-15) 


or 


M?(n)  -  M?(n_1) 


Mf 


(n) 


(B-1G) 


Typically  for  values  of  tj,  f2  ~  10  .  the  iteration  stabilizes  for  5  <  n  <  10. 


Difficulties  encountered  in  the  course  of  this  iterative  procedure  have 
definite  physical  significance  that  can  be  traced  to; 


a)  Either  of  the  characteristics  being  oriented  such  that  the  quantity; 

0  +  (-l)la  «  0  (B-17) 

L  Ji3 

or,  in  other  words,  a  characteristic  is  horizontal  in  the  flow  field 
for  which 

Q.„  -  «  •  (B-18) 

b)  Compressions  developed  in  the  flow  field  due  to  wave  coalescence 
such  that 

M?(n)  ^  1  .  (B-19) 

In  the  first  case,  the  quantity 

%  («3  -  Ri)  (B-20> 

must  be  reconsidered  when 

[o  +  (-l)loj  —  0  .  (B-21) 
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The  incremental  length,  M  ,  along  either  characteristic  is  given  by 


A £  - 


■  (S3-Ri) 


sin  0  -i 

(-l/a 

Ji3 

(B-22) 


and  as 


(-1) 


Hence,  the  quantity  in  the 


characteristics  equations  (B-ll)  and  B-12)  must  then  be  replaced  by  the  limit¬ 
ing  value  of  equation  (B-20),  i.e., 


sin  0  sin  cv 


(B-23) 


when 


0  +  (-1)  cv 


J  i3 


(B-24) 


The  second  problem  area  usually  results  from  "foldback"  of  the  charac¬ 
teristics  network  as  a  result  of  the  coalescence  of  the  same  family  waves-  To 
treat  this  problem,  provisions  must  be  made  in  the  overall  flow-field  claeulation 
sequence  so  that  the  shock  formed  in  the  flow’  field  as  a  result  of  this  coalescence 
can  be  treated  and  thus,  "foldbaek"  of  the  characteristics  network  avoided. 


3.  Axis  Points  (APS) 

If  any  one  of  the  points  (1,  2,  or  3)  is  located  on  the  axis  of  sym¬ 
metry,  then, 


R.,  0.  =  0  j  =  1,  2,  or  3  (B-25) 

and  the  calculation  procedure  of  section  2  is  modified  accordingly  for  the  axis 
point  calculations.  No  particular  problems  are  encountered  whenthe.se  condi¬ 
tions  are  imposed  on  the  field-point  calculation  sequence  since  the  term  [  R_1] 
only  appears  in  the  characteristic  equations  as  an  average  value  in  the 
coefficients  Q  . 

1  tJ 

The  two  axis  point  calculations  encountered  in  the  primary  flow-field 
analysis  are  shown  in  Figure  B-3. 


is 


FIGURE  B-3.  AXIS  POINTS 


For  the  one  case,  the  unknown  point  (3)  is  on  the  axis,  and  consequently, 

Ha,  03  =  0  .  (B-26) 

Then  the  remaining  values  [X,  M*]3  are  found  from: 


V  =  Yr 

(B-27) 

3  2  tan  (6  + 

a)  23 

M3  =  M3  -  P230  2  “ 

Q23P2  . 

'B-28) 

To  start  the  calculation  sequence,  it  is  assumed  that 

M3::  =  M*  . 

Then  a  first  approximation  [X,  M*l3  ^  can  be  found  from  equations  (B-27)  and 
(B-28).  Using  the  successive  approximation  technique,  the  calculations  are 
continued  until: 


(Mjf)  (Q)  -  (M3")(n"1) 


(B-29) 


Thus  the  conditions  at  point  (3)  are  determined  as:  [X,  0,  M*.  0] 


(n) 


In  the  other  case,  the  known  point  (2)  is  on  the  axis  where: 


R2,  0  2  —  0  • 

The  calculation  sequence  for  determining  the  values  at  point  (3)  are  the  same 
as  outlined  for  the  field -point  analysis  (section  2) . 


4.  Boundary  Points 


For  the  base-flow  analysis,  only  two  types  of  boundary-point  calcu¬ 
lations  occur,  viz. ,  the  constant -pres sure  condition  for  the  separated  flow 
region  and  the  solid  boundary  condition  when  an  afterbody  precedes  the  external 
stream's  separation  point. 


a.  Constant  Pressure  Boundary  Points  (CPBS) 

Along  the  boundary  (Figure  B-4)  the  condition  of  constant 
pressure  is  expressed  by: 

M*  ~  M|  =  M*  =  constant,  (B-30) 


where  M*  is  found  from  the  isentropic  flow  function 


P  =  constant. 
13 


(B— 31) 


R 
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A  first  approximation  to  the  spatial  location  of  the  unknown  point  (3)  is  found 
from  the  streamline  condition,  between  points  (2)  to  (3) ,  where 


(R3  -  R2)  =  (X3  -  X^)ianOa3 

and  along  the  I-characteristic,  between  points  (1)  to  (3),  where 
(R3  -  Ri)  =  (X3  -  Xi)tan  ( 0  -  a) i3  . 

Then  X3  is  found  from: 

x  _  [Ri  ~  R2  +  ^2  tan  6ri  -  Xt  tan  ( 0  -  n)  13] 

3  [tan  6  23  -  tan  ( 9  -  or)  13] 

The  local  flow  direction  is  ti  ~n  given  by: 


(B-32) 


(B-33) 


(B-34) 


0  a  -  0 1 


[  ( -  Mf)  -  Q13(R3  -  Ri)] 


t'is 


(B-35) 


The  calculation  sequence  is  initialized  by  assuming  that 


03  **  \  (01  +  0 2) 


The  corresponding  values  of  [x3^  ,  R3^ 


*(1)] 


are  then  found  from 


- - * - - - {_  o  >  — 0  »  o  j - - - 

equations  (B-30) ,  (B-34) ,  (B-32),  and  (B~35) .  By  use  of  these  equations  and 
the  current  approximations  to  the  variables  at  point  (3)  to  evaluate  the  average- 
value  coefficients,  the  estimates  for  [X,R,0]3  can  be  improved  by  successive 
approximations  until 


J  ( 


/’ 


(n) 


<  t 


Thus  the  variables  at  point  (3)  are  determined  as  [X,  Rs0]3in^  and  j^M3  =  iVL*J  . 


b.  Solid-Boundary  Points  (BTBPS) 


Along  a  solid  boundary  (Figure  B-5)  the  flow  must  be  tangent 
to  the  surface,  and  therefoi'e  at  any  point  (3)  on  the  surface,  the  local  flow 
direction  is  given  by: 


03  = 


(B-36) 


t 


Qj3  is  evaluated  accordingto section  2,  equations  (B-14)  or 


(B-23) . 
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FIGURE  13-5.  SOLID-BOUNDARY  POINT 


The  solid-boundary  profile  is  assumed  to  be  described  in  the  meridional  plane 
by  a  functional  relationship  of  the  form 


\  -  Rs  (X  )  •  (B-37) 


However,  before  the  streamline  condition  specified  by  equation  (B-36)  can  be 
applied,  an  initial  estimate  of  the  location  of  the  unknown  point  (3)  on  the  solid¬ 
boundary  must  be  made.  By  assuming  that  the  location  and  flow  properties  are 
known  at  a  point  (1)  in  the  adjacent  flow  field  and  a  point  (2)  on  the  solid¬ 
boundary,  initial  estimates  of  the  flow  properties  at  point  (3)  are 


O3  1  + 

M?  *  \  <M|  + 


9  2) 

(B-38) 

Mj )  . 

(B— 39) 

The  approximate  location  of  point.  (3)  onthe  solid  boundary  is  then  found  by 
solving  simultaneously  the  physical  I-characteristic  relationship  between 
points  (T)  and  (3)  and  the  equation  specifying  the  boundary  profile.  The  rela¬ 
tionships  to  be  solved  are: 
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(Rs  -  Rt)  =  (X3  -  Xi)tan  (0  a)  13 


(B-40) 


and 


Es  -  Ks  (Xs)  (B-37) 

when  R  =  R  and  X  -  X_.  Carrying  out  this  simultaneous  solution  and  using 

S  J  S  o 

equation  (B-36)  results  in  an  approximate  set  of  values  for 


[X,R,6]3(1) 


By  use  of  these  data  along  with  thd  hodograph  I-characteristic  relationship 
between  points  (1)  and  (3) ,  the  local  value  of  Mj  can  be  estimated  from 

M3  *  Mf  -  Pi3(63  -6  i)  +  PROD13  (B-4.1) 

where  for  1(0  -  u)13|  >  0  (cf.  Field  Point  discussion) , 


PROD  13  =  Q13(R3  -  Rt)  , 


(B-42) 


or  for  I  (0  -  o)  13 1  «  0 


PROD13  =  P14 


sin  0  sin  a 
R 


(X3  -  Xi) 


(B-43) 


Jl3 


Thus,  the  location  and  flow  properties  at  point  (3)  for  this  first  approximation 
are  determined  as: 


LX1R,M-,0]3(1)  . 

th 

For  the  next  and  all  successive  approximations,  the  (n  -  1)  estimates 

of  [M*,6]3  11  ^  are  used  with  equations  (B-40),  (B-37),  and  (B-41)  through 

(B-43)  to  determine  the  n^  estimates  of  [X,R,  M*,0  ]3^ .  This  procedure  is 
repeated  until  the  normalized  difference  in  Mg'  satisfies  the  convergence 
criterion: 


(B-45) 


t 


For  e  =  10-4,  equation  (B-45)  is  typically  satisfied  for  5  <  n  <  10. 
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Thus  when  equation  (B-45)  is  satisfied,  the  final  values  of 


[X,  R,  M*,0]3 


(n) 


are  determined. 


5.  Flow- Field  Analysis  (AC PBS) 

The  unit  processes  described  above  are  organized  into  a  sequential 
program  which  can  be  used  to  calculate  the  flow  field  subject  to  the  imposed 
boundary  conditions.  This  organization  is  principally  one  of  "bookkeeping’'  and 
is  normally  not  difficult.  The  basic  sequence  can  take  any  of  several  different 
forms;  the  preference  here  has  been  to  calculate  along  Family  I  characteristics 
toward  the  boundary. 


a.  Calculation  Sequence 


The  boundary  conditions  imposed  on  the  flow-field  calculation 
(  Figures  B-6  through  B-9)  are: 

1)  The  flow  variables  are  specified  along  the  initial  internal 
(nozzle)  flow  or  external  (free  stream)  flow  characteristic. 

2)  The  conditions  along  the  boundary  are  specified  for  the 
particular  type  of  boundary  condition  that  is  being  considered. 

These  data  along  with  the  surface  geometry,  if  applicable,  are  sufficient  to  deter¬ 
mine  the  flow  field. 

At  the  nozzle  corner  ^X^,  R^a  centered  expansion  (or  compression), 

can  occur  as  a  result,  of  the  need  for  the  internal  flow  to  expand  (or  compress) 
at  the  nozzle  exit  to  satisfy  the  imposed  constant-pressure  boundary  condition. 

A  similar  situation  can  also  occur  for  the  external  (free  stream)  flow  as  a  result 
of  expansion  (or  compression)  corners  and  a  solid-boundary  preceding  the 
separation  point  of  the  external  stream  as  well  as  the  requirement  that  the  exter¬ 
nal  flow  must  adjust  to  the  pressure  level  maintained  downstream  of  the  separa¬ 
tion  point.  The  general  flow- field  calculation  sequence  selected  here  for  Uher 
the  internal  or  external  flow-field  proceeds  from  the  initial  Family  II  charac¬ 
teristic  (nozzle  or  afterbody)  along  I-charaeteristics  to  the  boundary. 

The  calculation  sequence  for  internal  flow  and  external  flow  are  illustrated 
in  Figures  B-C  and  B-7,  respectively. 
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INITIAL  1 1- CHARACTERISTIC 


FIGURE  B-6.  FLOW- FIELD  ANALYSIS 


b.  Wave  Coalescence 

The  Family  II  characteristics  from  the  internal  constant- 
pressure  boundary  tend  to  steepen  as  the  calculations  proceed  in  the  downstream 
direction.  These  characteristics  eventually  coalesce  and  form  a  shock  wave 
within  the  flow  field.  This  condition  is  detected  by  the  crossing  of  waves  of  the 
same  family  thus  giving  rise  to  the  "foldback"  of  the  characteristics  network. 
Although  flow-field  calculations  where  "foldback"  occurs  still  yield  results 
which  are  in  reasonable  agreement  with  experiment  [7] ,  the  flow-field  calcula¬ 
tions  must  invariably  be  terminated  as  a  result  of  errors  directly  attributable 
to  the  unrealistic  characteristic  network  that  develops. 

An  exact  treatment  of  this  problem  has  been  given  |8] ,  as  well  as  an 
approximate  treatment  [9].  However,  since  the  calculated  boundary  is 
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FOLD  BACK 


FIGURE  B-9.  MODIFICATION  OF  THE  CALCULATION  SEQUENCE 

■PA13  n-'TVC' A 'T'Tyrv*  xita  171?  nAA  r 
a  val  A  XUjA  A  liXU  KlfnVii  UOA  A-lL^OV  Ai 
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relatively  insensitive  to  the  method  ir  en  10  treat  the  coalescence  problem,  a 
simplified  approach  will  be  described  here  which  is  in  excellent  agreement  with 
previous  reports  [  s,  y] . 

In  reference  to  Figure  B-8,  wave  coalescence  has  occurred  within  the 
flow  field.  Allowing  a  single  "foldback"  at  this  point,  the  conditions  on  the 
"upstream"  and  "downstream"  sides  of  the  coalescence  point  are  determined  by 
linear  interpolation  between  the  points  (2)  and  (3)  and  the  points  (2)'  and  (3)', 
respectively.  The  "foldback"  points  (3)  and  (3)'  are  then  dropped  and  the  flow- 
field  calculations  are  continued  by  using  the  flow  variables  determined  at  the 
wave -coales cent  point.  Figure  B-9. 

In  actuality,  the  oblique  shock  wave  formed  in  the  flow  field  due  to  the 
wave  coalescence  propagates  downstream  where  it  becomes  curved  as  a  result 
of  the  continuous  interaction  between  the  shock  and  the  waves  in  the  flow  field. 

As  a  consequence,  the  flow'  is  rotational  downstream  of  the  internal  shock  wave. 

In  the  approximation  described  above,  the  flow  is  assumed  to  remain  irrotational. 
This  assumption  yields,  in  most  cases,  results  which  are  acceptable  and  consis¬ 
tent  with  the  overall  analysis. 


6.  Initial  Nozzle  Characteristic 


The  nozzle  geometries  are  restricted  to  those  configurations  which 
produce  sonic,  uniform,  or  conical  supersonic  flow.  The  objective  is  to  deter¬ 
mine  the  flow  variables  along  the  initial  nozzle  characteristic  for  each 
configuration. 


a.  Sonic  Nozzle 


The  sonic  nozzle  can  be  treated  approximately  as  a  nozzle  which 
produces  unifomi  flow  at  the  nozzle  exit  that  is  slightly  supersonic,  e.g. , 


b.  Uniform  Supersonic  Nozzle  (UFLOC) 

For  uniform  supersonic  flow  at  the  nozzle  exit  (Figure  B-10) 
the  initial  I-charaeteristic  is  straight  and  the  flow  variables  are  known  as 
[Mu. '!  II  ^  a*on&  this  characteristic. 

The  case  where  an  initial  compression  must  exist  at  the  nozzle  exit  to 
satisfy  the  imposed  boundary  condition  will  be  considered,  in  an  approximate 
way,  in  the  subsection  to  follow,  ”A  Compression  at  the  Nozzle  Exit. " 
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c.  Conical  Supersonic  Nozzle  (CNFLOC) 

The  flow  in  an  ideal  conical  nozzle  (Figure  B-ll)  can  be 
specified  as  being  at  a  uniform  Mach  number  ^or  along  the  zone  of  the 

spherical  sector  that  coincides  with  the  nozzle.  The  initial  I-characteristic  can 
be  determined  exactly  [10],  or  numerically  from  the  flow  conditions  specified 
on  the  non-characteristic  spherical  surface. 


The  numerical  approach  while  yielding  values  for  the  initial 
I-characteristic  that  are  in  excellent  agreement  with  previous  exact  theoretical 
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results  1 10  i  in  a  tin-  advantage  of  being  easily  extended  to  treat  approximately 
the  problem  of  an  initial  compression  at  the  nozzle  exit  (for  both  the  uniform 
and  conical  nozzles)  . 


i'he  flow  conditions  are  specified  along  the  non-characteristic  curve 
(  Figure  11—  1 2)  as  being  at  a  uniform  value  of  M*  and  that  the  velocity  vector  is 


always  perpendicular  to  this  curve. 

where  o  s  o  ^  <j 
JNC  11 


X.R.M-yd 


Any  point  on  this  curve  is  defined  by 
and 


X  =  X.  +  R  fcos  0  -  cos  0 ,  1 
II  cone  [  llj 


R  =  R  sin  0  , 
cone 


where 


R  =  Rlr /sin  0  . 

cone  11/  II 


(B-46) 

(B-47) 


(B-48) 


FIGURE  13-12.  CHARACTERISTICS  NETWORK  FOR  NUMERICAL 
CALCULATION  OF  THE  INITIAL  II-CHARACTERISTICS 
LOR  A  CONICAL  NOZZLE 


The  non-characteristic  curve  is  then  subdivided  and  a  characteristics 
network  (  Figure  B-12)  can  be  used  to  determine  numerically’  the  flow  variables 
and  corresponding  locations  on  the  initial  conical  nozzle  characteristic. 

This  general  calculation  sequence  is  useful  since  the  calculations  are 
made  from  a  non-characteristic  curve  to  the  corresponding  characteristic  curve. 
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d.  A  Compression  at  the  Nozzle  Exit 


If  the  boundary  conditions  are  such  that  the  nozzle  discharges 
into  a  region  at  a  pressure  greater  than  the  supersonic  design  pressure,  an 
oblique  shock  wave  is  generated  in  the  primary  flow  field.  If  this  compression 
is  assumed  to  be  relatively  weak,  the  oblique  shock  can  be  treated  approxi¬ 
mately  as  a  reversible  compression.  For  either  the  uniform  or  conical  nozzle, 
the  flow  variables  on  the  noncharacteristic  curve  and  the  imposed  boundary 
conditions  are  used  to  establish  a  single  reversible  compression  wave  at  the 
nozzle  exit  location  R^  which  satisfies  the  boundary  conditions.  With 

these  data,  the  remainder  6f  the  initial  I-characteristic  can  be  established 
using  the  calculation  sequence  from  the  non-characteristic  curve  to  the  corres¬ 
ponding  characteristic  curve  (Figure  B-13), 


FIGURE  B-13.  APPROXIMATE  ANALYSIS  OF  A  COMPRESSION 
AT  THE  NOZZLE  EXIT 


7.  Afterbody  Analysis  (ABTS) 

The  afterbody  analysis  utilizes  the  unit  processes  described  in 
section  4;  this  analysis  must,  however,  be  consistent  with  the  flow-field  analysis 
of  section  5a  (Figure  B-7) .  In  addition  to  determining  the  detailed  flow  condi¬ 
tions  over  the  afterbody,  the  final  Family  II-characteristic  originating  at  the 
terminus  of  the  afterbody  must  be  determined  so  that  the  same  flow-field  calcu¬ 
lation  algorithm  can  be  used  for  calculating  both  the  external  and  internal  flow 
fields  and  constant-pressure  boundaries. 


31 


Initial  and  Final  Family  II  Characteristics 

The  approach  flow  to  the  afterbody  is  currently  restricted  to 
uniform  supersonic  flow.  For  this  case,  the  initial  II-characteristic  is  a  straight 
line  along  which  the  flow  properties  are  constant  and  known  at  each  arbitrarily 
selected  subdivision  point;  this  situation  is  shown  in  Figure  B-14. 


FIGURE  B-14.  INITIAL  II-CHARACTERISTICS  FOR  UNIFORM 
EXTERNAL  SUPERSONIC  FLOW 


When  an  afterbody  is  present,  the  points  along  the  final  II-characteristic 
originating  at  the  afterbody  terminus  are  determined  as  input  data  for  the  sub¬ 
sequent  constant-pressure  boundary  calculations.  After  the  I-characteristi  ■ 
originating  at  the  initial  II-characteristic  and  passing  through  the  terminus  of 
the  body  is  found  by  iteration,  the  remainder  of  the  points  along  the  final 
1 1- characteristic  are  found  by  continuing  to  subdivide  the  initial  II-characteristic 
and  continuing  the  method  of  characteristics  calculations  to  the  final 
II-characteristic  (Figure  B-15). 

b.  Specific  Afterbody  Profiles  (BTCNST) 

As  discussed  in  section  4b  ("Solid-Boundary  Points") ,  a 
simultaneous  solution  must  be  made  between  the  equations  describing  the 

:i2 


FIGURE  B-15.  DETERMINATION  OF  THE  FINAL  II-CHARACTERISTIC 
THROUGH  THE  AFTERBODY  TERMINUS 


afterbody  profile  and  the  physical  I-characteristic;  this  solution  determines 
points  in.  the  characteristic  network  which  are  located  on  the  afterbody  surface. 

The  purpose  of  this  section  is  to  summarize  the  resulting  expressions 
for  three  afterbody  shapes  —  the  ogive  (circular  segment) ,  the  parabola,  and 
the  cone.  A  typical  afterbody  configuration  and  the  associated  notation  are 
given  in  Figure  B-16;  for  each  shape,  the  vii-rnody  is  assumed  to  be  completely 
specified  by  the  values  of 

RBT1’  ^BTlj  aUd  [XBT2*  RBT2]  ‘ 

The  equations  specifying  the  afterbody's  meridional  profile  are  sum¬ 
marized  in  Table  B-I.  The  constants  [CijCjjCg]  in  the  profile  equations  aie 
determined  for  each  afterbody  based  on  the  values  specified  at  points  (BT1) 
and  (BT2).  The  resulting  expressions  for  these  constants  are  summarized  in 
Table  B-H. 
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FIGURE  B-16.  GENERAL  AFTERBODY  NOTATION 


TAB LE  B-I.  EQUATIONS  OF  AFTERBODY  MERIDIONAL- 
PLANE  PROFILES 


Shape 

R  =  R  (x  \ 

s  s  \  s/ 

Ogive 

Parabolic 

Conical 

■■  - - —  -  ■  _ _ 

'.-'i'  [c3-  (xs  - c 2)^ 
Bs  =  c3  +  CA  +  °1XS 

The  simultaneous  solution  of  equations  (B-37)  and  (B-40)  has  been 
carried  out  for  the  profiles  specified  in  Table  B-l;  the  r-esults  of  these  solutions 
are  summarized  in  Table  B— III. 


TABLE  B-II.  CONSTANTS  IN  AFTERBODY  PROFILE  EQUATIONS 


i 


TABLE  13— III.  EXPRESSIONS  FOR  [X3l  R3>  fl3)  FROM  THE  SIMULTANEOUS 
SOLUTION  OF  EQUATIONS  (B-37)  AND  (B-40)  FOR  VARIOUS 
AFTERBODY  SHAPES 


— 

Shape 

Simultaneous  Solution  Values  for  [X3l  R3,  03] 

Ogive 

v  _  !"  B  -[B2  -  4A  C  ]  ^2  ) 

2A 

A  =  1  +  tan2  (9  -  a)  13 

B  =  2(R!  -  C[)tan  { 0  -  Qf)13  -  2C2  -  2Xj  tan2  (d  -  a)  13 

C  =  C22  -  C3  +  I  (R3  -  Ct)  -  Xt  tan  (6  -  a)l3]2 

R3  =  Rj  +  (X3  -  X^tan  (9  -  a)iS 

Parabolic 

„  -  B  -[B2  -  4AC]/2  } 

2  A 

A  =  Ct 

B  =  C2  -  tan  (0  -  a)13 

C  =  C3  -  Rj  +  Xt  tan  (0  -  o)13 

R3  =  Rj  +  (X3  -  XJ  tan  (0  -  a)  iS 

03  =  tan"1  [C2  +  2CjX3) 

Conical 

v  _  [C[  -  Rj  *•  C2C3  +  Xi  tan  (9  -  o)i3] 

'  3  [tan  (0  -  aO  13  -  C2] 

R3  —  Rj  +  ( X3  —  ( 0  —  cy)  13 

03  =  tan-1  (Cj)  =  constant 


Appendix  C 

MIXING  ANALYSIS:  TWO-DIMENSIONAL 
CONSTANT-PRESSURE  TURBULENT  MIXING 


The  general  references  for  Appendixes  C  and  D,  although  not  cited,  are 
several  by  Korst  et  al.  [11-16]. 


1.  Two-Dimensional  Constant- Pressure  Mixing  Region 

The  two-dimensional  mixing  region  control  volume  is  shown  in 

Figure  C-l. 


For  constant-pressure  mixing  with: 

a)  uniform  flow  at  x  ~  0 

b)  tully  developed  flow  at  x,  the  velocity  profile  within  the  mixing  zone 
is  given  by 


y  y 


0  =  7J-=|  [1  +  erf  7)]  , 
a 


where 


in  which 

(x,y)  refer  to  the  intrinsic 
coordinate  system 

o  =  similarity  parameter 

a  =  12  +  2.  76  M  . 

a 


FIGURE  C-l.  TWO-DIMENSIONAL  MIXING 
REGION  CONTROL  VOLUME 


2.  Determination  of  y 

m 


The  intrinsic  coordinate  system  is  located  relative  to  the  reference 
coordinates  by  applying  the  momentum  equation  in  the  x-direction  (per  unit 
width)  to  the  control  volume  of  Figure  C-l.  The  result  is: 
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Now  at  V  anu  ’W 


|1  -  4>l  <  <a*  | 

IT  -  T  I  <  e  " 
oa  o  a 


for  rj  —  tj 


ra 


for  „  *  VRb 


I  4>\  <  V 

Itb  -  V  <  S' 

Therefore,  equation  (C-3)  can  be  rearranged  as 

Sa  (7,m+!7Ra) 

-V+/,  „  4^=0 

(V%b)  a  ^  a 


'Ra 


Thus 


'Ra 


+  ^  —  <f>2dn 


\b 


<Jfd  f)  -*-77  -  0 


m 


and  consequently, 


'Ra 


ra 


-  '  /  p-  ^ 

"Rb  a 


3.  Identification  of  j-streamline 

Application  of  the  continuity  equation  to  the  jet  flow  yields 


Ra 


p  U  Y  .  +  f  P udY  =  0  . 

a  a  RA  J 


( C— 4) 


(C-5) 


(C-6) 


( C— 7) 


(c~s: 


(C-9 


In  non-diinensional  form,  equation  (C-9)  becomes. 


(C-10) 


or 


y% 


"W  J 


P  u  , 

~  u*  d" 

a  a 


Ra 


-’Wf  J  7  4,d"+  J 


>Jj  a 


yv 

— -  <pdrj  =  0  , 
P„ 

7L-  2 

'Ra 


(C-ll) 


T)Ra 

J 

ri5 


P 


<Sdrj  = 


( C-12) 


anti  combining  equations  (C-8>  and  (C-12),  the  integral  equation 
71.  is 


7W 

J 

v, 


<f>dri 


7)Ra 

/ 


be  solved  for 


(C-13) 


4.  Energy  Considerations 

lly  use  of  Crocco's  Integral  Solution  for  an  apparent  unity  turbulent 
Prandtl  number,  the  stagnation  temperature  distribution  in  the  mixing  region 
is  given  by 


T 

— -  =  A  + 
Toa 


oa 


( C— 14) 


( C— 15) 
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or 


where 


A  '  Ab  +  (‘  -  ab)*  • 


A 


B 


a*  Energy  Transferred  Across  j-Streamline 
Shear  work  and  conduction  per  unit  width  is 


,(**♦»») 

J  PC 


p 


In  dimensionless  form, 


p  -j  (W%) 

xpUCT  “  J  (1 

a  a  p  oa  7j. 


since 


(1  -  A)  a  0  for  n  <  rj  < 


("a,  *  ”n,) 


(C-1C) 


( C — 1 7 ) 


(c-ia) 


(c-ia) 


(C-20) 


equation  (C-20)  ear,  be  simplified  as-. 
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b  !  X 


f? 

ir  j 
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v  r  r  t  J  {1  -  A>  T  ^d,< 

X  p  l  t  i  J  p 

a  a  p  oa  //^  a 


'ha  "ila 
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Hy  use  of  equation  (C-lfl), 
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Also  from  equation  (C-13), 
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Therefore, 
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becomeltoUingthe  Simplification-  ^  dimensionless 


energy  transfer  rate 


£ _ 1 

x  p  C  U 


T  =-[ 


E  P  a  °a  r4tb  a 


B  *  n 
?Rb  a 


^77 


’  {  R1  ‘  Ab)*]f  ***  ■ 

rib  J  a 


(C-  21) 


From  equation  (C-17) 

A  =  AB  +  (X  *  ab)  0  or 
hence,  equation  (C-21)  can  be  expressed  as: 


(i-ab)»].(a-ab)  . 


x  p  u  C  A 

a  a  p  oa 


T  J  (A  ~  ab)  ~  <pdr) 

oa  n  '  pa 


(C-22) 


'itb 
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b. 


Energy  Convocted  with  Maas  Between  (j )  'and  (d)  Streamlines 
(with  Reference  to  7,e ro  Temperature)  per  Unit  Width 


The  energy  eonvected  ia: 


e  ,  -  I  pC  uT  dy 
a  J  p  o  J 


(C-23) 


nr  in  dimensionless  form 


,  % 
d  =  /  f  A^dr, 


x  p  C  U  T  J  p 
a  p  a  oa  n.  a 
J 


and  finally  expanding  the  integral 


?,d  7?Rb  , 

—  - =  j  --  A 0d r?  +  j  -  A0dp  . 

x  p  C  IT  p  J  p 

a  p  a  oa  ^  a  V.  a 


(C-24) 


11  a  ;  -  n  *  i  then 


<j\  -  -d7] 

and  the  last  integral  in  equation  (C-241  can  be  exp  re  a  a  ed  aa: 


,JRb  ^ 

J  A0d„  =  -  /  ^A0dX 

’j  1  %b  a 


since 


4>(rj)  —  0 


(”»>  +  ")) 


-  \ 


>  0  for  n,  <  \  <  77. 

Rb  j 


4-1 


and 


A  >  (>  , 


Therefore,  the  dimensionless  convected  energy  transfer  rate  is  given  by: 


x  p  C  U  T 


~  J  7T  "  I 


A  (/'dp  . 


(C-25) 


a  p  a  oa  a  ^  a 


The  total  rate  of  energy  transfer  per  unit  width  to  the  wake  is  found  by 
combining  equations  (C-21)  and  (C-25).  The  result  in  dimensionless  form  is 


x  p  C  U  T  xp  C  U  T 
a  p  a  oa  a  p  a  oa 


(eJ  +  *«) 


(C-26) 


x  p  C  U  T 
a  p  a  oa 


V. 

-A  J  J—  <pdri  +  f  — 

B  J  l>  J  p 

\ib  a  %b  a 


A(pdi) 


(C-27) 


5.  Mass  Transfer 


Mass  rate  convecied  between  (j)  and  (d)  streamlines  per  unit  with 


width  is 


=  J  P  u  dy 


y. 

j 


(C-28) 
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or  in  dimension  loss  form 


xp  u 
a  a 


d  ] 

J  J~  W*  -  I  J~  ^  ■ 

?iRb  a  %b  a 


(C-29) 


Equations  (('-  12),  (C-27),  and  (C-29)  define  results  for  mixing  analysis  after 
4>and  —  are  determined. 


(j.  Constant- Pressure  Mixing 

For  an  ideal  gas,  the  equation  of  state  is 

P  =  piRT  ; 

then  for  P  =  constant, 


p  T  =  pT 
a  a 


(C-30) 


or 
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From  the  energy  equation 
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Defining  the  free  stream  Crocco  number  squared  as 


U  2 

p  2  — _ ^ 

a  “  2C  T 

p  oa 
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—  -  0 

a 

T  +  C  202T  =  T  ,  (C-33) 

a  oa  o 


at  "a,  ” 


Therefore, 


at  any  </>, 


T  +  C2  T  =  T 
a  a  oa  oa 


(A  -  C.V) 


Therefore,  from  equations  (C-31),  (C-34),  and  (C-35) , 


_P 

P 

a 


(034) 


(C-35) 


(036) 


(C-37) 
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Defined  Integral  Notation 


The  integral  notation  is  defined  by  the  following  equations: 


1 1  (  ?7> 


■HL) 


(A  -  C^) 


(C-38) 


i2  ( »;) 


J  ( A  -  C  W)  d ''  '  12  (n’  ’  Ca2) 

7^tib  V  a  ) 


(C-39) 


where 


( C-40) 


s,  liquation  Summary 

The  equation  to  be  solved  for  u.  is: 

''itu  /'j  >fRa  ?)Ra 

I  —  <f)dr;  =  -  j  —  (pdrj  +  I  —  <jt>d 77  =  I  —  <f>2dr)  (C-13) 

J  0  J  p  *'  p  J  p 

"j  a  "lib  a  %b  a  %b  a 

or  in  terms  of  the  defined  integrals,  equations  (C-38)  and  (C-39) 

‘ti)- ['■(%>)  -  H’^)]  • 

The  relationship  between  I3  and  11#  I2  is 


(C-41) 


or  in  terras  of  the  defined  integral  notation, 


k  -  abi,  +  (i  -  a^i,  . 

From  equation  (C— 27) ,  the  energy  transfer  rate  becomes 

ere _  r 

ph^aV^oa  ‘  IK’x)  -  V'(”j) 

From  equation  ( C— 29) ,  the  entrained  mass  rate  is 

e  -4'aI1’  W" 11  w.  ■ 


(C-42) 


(C-43) 


9,  d-St  ream  line  Velocity  Ratio 


For  the  d-streamline,  the  recompression  criterion  yields  a  value  of 

the  d-streamline  Crocco  number,  C  The  value  of  C  ,,  the  adjacent  free  stream 

d  d 

data,  and  the  state  of  the  quiescent  fluid  can  then  be  used  to  determine  the 
d-streamline  velocity  ratio  <^.  This  is  accomplished  in  the  following  way. 


The  definition  of  the  velocity  ratio  is  considered: 

S  “d/Ua 

or  rearranging  in  terms  of  the  corresponding  Crocco  numbers, 
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( C— 44) 


(C-45) 


The  stagnation  temperature  ratio  can  be  expressed  from  equation  (C-16)  as 


where 
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Ther«  i'ore  by  eombinintt  equations  (C-4f>)  and  (C-4G),  the  expression  to  be 

solved  tor  />  is 
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-  ab)  ‘V 
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where  the  ratio  of  the  Crocco  numbers  has  been  defined  as 


Appendix  D 


MIXiNG  ANALYSIS:  APPLICATION  OF  THE 
TWO-DIMENSIONAL  MIXING  ANALYSIS  TO  THE 
TWO-STREAM  AXISYMMETRIC  BASE-PRESSURE  PROBLEM 


The  approximate  two-stream  flow  model  is  shown  in  Figure  D-l. 


(IE)  PBE 


FIGURE  D-l.  APPROXIMATE  TWO-STREAM  FLOW  MODEL 

I 


1  *  Two-Dimensional  Approximation  of  the  Mass  and  Energy 
Transfer  Rates 


Superposition  of  the  two-dimensional  mixing  region  on  the  "corre¬ 
sponding"  inviscid  flow  field  at  the  impingement  point  (Appendix  C)  of 
Figure  D-l  is  shown. 


a.  External  Stream 


equation 


The  approximate  entrained  mass  tiow  rate  is  determined  from 
(C-43)  as 


C  =  2ttR. 


i.41  -  cbk)-4',ueuiJe[I|(V:)  -  ‘‘(V:)]  '  tU_1> 
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Ih<-  approximate*  trunsfor  rate  is  determined  from  equation  (0-42>  as: 


i  mp 


(’  ‘  1  lu.)  i^7  'hi:1  m^pidoK  l3('(ii:)  ~  'W:*1 ^ji:) 


win* iv  A  T  1' 

ME  M  oE 


h.  Internal  Stream 

The  ent  rai  ihh]  muss  flow  rate  is  (equation  C— 13) 

S 

1 


G.  -  ‘J-H 


imp 


0  "  1  Ml)  <x,  '’lJl1  Ml 


The  energy  transfer  rate  is  (equation  C-42) 


(D-2) 


( D-:i) 


S 

■  -r"mP ('  -  ’  V>(v)l  ■,D-W 


where  A  -  T,,  T  . 
Ml  11  i)l 


2 .  Gonservation  of  Mass 

For  the  base-region  euntrol  volume,  conservation  of  mass  requires 

that 


G  >  G,.  -x  G,  -  0  . 
o  E  1 


(D-5) 


Tile  term  li  ^  is  a  mass  flow  rale  entering  the  base  region  at  a  ncitligible  velocity, 

i.c.,  "mass  hlci’il.  "  My  introduction,  as  a  reference,  of  the  internal  stream's 
f nozzle  lion.)  mass  flow  rate,  G  ,  equation  (I)-'i)  can  be  written  in  dimension- 
less  lorin  as 


whore  the  "bleed  ratios"  are  defined  by 


U.i  :  t:o/Sl 

5  :  -  H  «,)/uki  ' 


(D-7) 


(D-x) 


li.  Conservation  of  Energy 


For  the  base  region  control  volume,  conservation  of  energy  requires 

that 

E  +  /E  -t  E 
o  V  E 

By  use  of  the  total  enthalpy  of  the  internal  stream  (nozzle  flow)  as  a  reference, 
i.e.  , 


=  0 


(D-D) 


EN1  '  ^NI^pl^ol  ’ 

equation  (D-S)  can  be  written  in  dimensionless  form  as 

K  -  E  -  0  , 
o 


(0-10) 


(D-ll) 


where 


£„  S  V‘-NI  IU-12) 

E  :  -  (H  +  e,)/kni  • 


Thus  for  conservation  of  mass  and  energy  in  the  base  region,  equations  (D-U) 
and  (D-ll)  must  be  simultaneously  satisfied. 


4.  Nozzle  Mass  Flow  Kate,  ti 


Figure  U-2  presents  the  ideal  nozzles  and  their  notation. 


For  source  flow  through  the  eonieal  nozzle  (Figure  D-2a)  the  flow 
Mach  number  is  a  constant  over  the  nozzle's  exit  area  and  the  velocity  vector  is 
always  perpendicular  to  this  area.  The  exit  flow  area  is  equal  to  the  area  of  the 
zone  of  a  sphere  of  radius  K  .  The  exit  area  is  found  as: 


who  re 


m  =  gas  constant  (lbf-tt/lbm-QU) 

A  =  exit  1'low  area,  equation  (D--14)  {  L2) 

P  =  exit  pressure  (lbi/L‘) 

T  =  stagnation  temperature  ('  ll) 

ol 

g  =  32.174  (Ibm-i't/lbi-sec2) 
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<J  =  mass  flow  rate  (lbm/sce) 
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then 


Vu 


The  function  F  ( M,  y)  can  be  expressed  in  terms  of  M*  as 


f<m*.y)  = 


2v 


v  <  i 

or  the  nozzle  mass  flow  rate  becomes 
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(D-20) 


(D-21) 


(D-22) 


(D-23) 


For  uniform  flow,  (3  =  0,  and  conical  flow  l>  >  0. 


5.  The  Quantity  (pU) 


u 


For  an  ideal  gas, 


(D-24) 


For  3R  in  units  of  ( Ibf-ft/lbm- ° R) ,  (pU)  is  given  by 


(pU) 


(D-32) 


where  in  equation  (D-32) 

p  (lb in/ ft3) 

U  (ft/sec) 

P  (Ibi/  L2) 

T  ( °R) 
o 

H  (lbf-ft/lbm-uR) 
gc  (Ibm-ft/lbf  -secz) 


(i.  Dimensionless  Bleed  Ratio,  B 


The  dimensionless  bleed  ratio  [equation  (D-8)  ]  is 


(0-8) 


By  substitution  of  equations  (D-l) ,  (D-3),  (D-23),  and  (D-32)  into  equation 

(D-8)  and  by  using  the  condition  that  P„  =  P_  T  =  P.  the  result  is 
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7.  DimenBioniess  Energy  Ratio,  E 


(D-33) 


The  dimensionless  energy  ratio  [equation  (D-13)  ]  is 

E  =  -  (ee  +  Ei)/eni  . 


(D-13) 
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I.V  HuioUtutnui  ill  imi;;  ,  (1>-  l) ,  (L)-10),  ( 0-33)  ,  ;uul  (1J-I52)  into 
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b.  (' well ieie tits:  :i|,  a2<  ;i;i 

To  simplify  the  notation  of  equations  (D-I3'i)  and  (D-34) ,  the  follow¬ 
ing  coeifieu'iits  have  been  delined.  They  are; 
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8iiii|ilii'iu(i  Forms  oi  (13,  K) 


From  equations  (D-G'.S),  (D-GGl,  and  (l)-:sf>)  ,  the  dimensionless 
nia-^s  bleed  ratio  becomes 
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and  from  equations  ( D-:i4) ,  (D-G5) ,  and  (D  -:t7),  the  dimensionless  energy 
transi’er  ratio  becomes 
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